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Abstract 

A microcanonical finite-size scaling ansatz is discussed. It exploits the existence 
of a well-defined transition point for systems of finite size in the microcanonical 
ensemble. The best data collapse obtained for small systems yields values for the 
critical exponents in good agreement with other approaches. The exact location of 
the infinite system critical point is not needed when extracting critical exponents 
from the microcanonical finite-size scaling theory. 
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In the canonical ensemble phase transitions appear exclusively in infinite sys- 
tems. At the critical point of a continuous phase transition denoted by T^ and 
he the Gibbs free energy as a function of the temperature T and of the field 
h conjugate to the order parameter develops a non- analytic behaviour. As a 
consequence some thermodynamic functions as e. g. the specific heat and the 
susceptibility exhibit power law singularities \T — Tc\~" and \T — Tc\~'^ in the 
vicinity of Tc and h = he fixed, commonly with non-classical critical exponents 
a and 7. In finite systems all these singularities are rounded in the canonical 
ensemble. The appearence of this rounding is taken into account by finite-size 
scaling theory originally proposed on phenomenological grounds [1,2,3]. In the 
asymptotic limit L ^ 00 and T —^ Tc the behaviour of finite-size quantities is 
governed by scaling functions. The scaling functions are basically determined 
by the ratio L/^{T) with ^(T) being the correlation length and L the linear 
extension of the system. Note that the validity of finite-size scaling relations 
has been proven in renormalization group theory [1] . Over the years finite-size 
scaling theory has been shown to be a valuable tool for extracting critical 
exponents from finite-size data, as obtained for example from numerical sim- 
ulations. 

In the microcanonical analysis of finite systems one considers the entropy 
Sn{E, M) = lnQi\f{E, M) as a function of the energy E and of the magneti- 
zation M [4]. The density of states of a system with A^ spins is denoted by 
Qn{E,M) and natural units with ks = I are chosen. The entropy surface 
exhibits a well-defined transition point at an energy Ec and magnetization Mc 
although in finite systems Sn{E, M) is everywhere perfectly analytic. The mi- 
crocanonical analysis also shows that typical features of symmetry breaking, 
as for example the abrupt onset of several order parameter branches when 
the transition point is crossed from above, are already encountered in small 
systems [5,6,7]. With regard to these intriguing effects, it is tempting to ask 
whether a direct analysis of the microcanonical entropy also allows the deter- 
mination of critical quantities from finite-size data. 

In this Letter we provide numerical evidence that it is in principle possible to 
extract critical exponents from the behaviour of finite microcanonical systems. 
Considering different classical spin models belonging to different universality 
classes we show that the critical exponents can be obtained from systems 
which are moderately small. The scaling ansatz to be discussed in the following 
does not need any a priori knowledge of the infinite system. Especially, the 
knowledge of the exact location of the critical point of the infinite system is 
not needed. 

In the following we discuss the nearest neighbour ferromagnetic Ising (I) model 
in three dimensions with the Hamiltonian 

nj = -jY,s,Sj (1) 



as well as a generalized Ising (GI) model with equivalent nearest and next- 
nearest neighbour interactions defined by the Hamiltonian 

Hgi = —J '^ SiSj — J ^ Si Sk- (2) 

Here, J > is the coupling constant and the spins Si can take on the val- 
ues d=l. The first sum in equation (2) extends over nearest neighbour bonds 
whereas the second sum is over bonds connecting next-nearest neighbours. On 
the simple cubic lattice, a spin has 6 nearest neighbours and 12 next-nearest 
neighbours which yields the ground state energy per spin Egi = —9 J for the 
GI model. On the same lattice the simple Ising model has the ground state 
energy ej = —3 J. The continuous phase transitions observed in both models 
belong to the same universality class. In addition, we investigate the three 
state Potts model in two dimensions whose Hamiltonian reads 

np = -jY.5s^,s,, (3) 



with the spins taking on the values Si = 1, 2, 3. 

In a microcanonical analysis of finite systems with N = L'^ spins {d being the 
number of space dimensions), the object of interest is the density of states 
Qn{E, M) as a function of the energy E and the magnetization M. The spon- 
taneous magnetization Msp^n{E) is defined to be the value of M where the 
entropy Sn{E,M) at a fixed value of the energy E has its maximum with 
respect to M : 

M,p,;v(^) : ^^ Sn{E,Msp,n{E)) = niax5^(E,M) . (4) 



This definition of the microcanonical spontaneous magnetization assures the 
equivalence of the canonical and microcanonical ensemble in the thermody- 
namic limit. At energies lower than a transition energy E^^n the entropy of 
the Ising universality class exhibits two maxima at M = ztM^pAr. When ap- 
proaching the finite-size transition point E^^n from below the order parameter 
Msp^N vanishes with a square root behaviour in the finite system [5,7]. Thus 
the transition energy Ec^n can be localized with high precision. This vanishing 
is also refiected in the divergence of the microcanonically defined finite-size 
zero- field susceptibility [5]. In Figure 1 the variation of the spontaneous mag- 
netization per spin, m^p tv = Mgp^j^/N, as a function of the specific energy, 
e = E/N, is shown for the two-dimensional nearest neighbour Ising model on 
an infinite and on a finite square lattice with 32^ spins. Both curves coincide 
at low energies. At higher energies, the square root behaviour close to the 
finite-size transition point can be observed. A similar behaviour is found for 



the entropy surface of finite two-dimensional three state Potts models [8]. Note 
that the discreteness of the physical quantities of discrete spin systems such 
as the Ising model has to be considered with some care. There the language 
used here refers to continuous functions that describe the discrete data most 
suitably. In the microcanonical analysis of continuous spin systems as e.g. the 
XY [9] or the Heisenberg model this concern does not exist. 
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Fig. 1. Microcanonically defined order parameter vs specific energy for 
two-dimensional Ising models defined on an infinite square lattice (full line) and 
on a finite square lattice containing 32 x 32 spins (symbols). The finite-size sponta- 
neous magnetization vanishes at a well-defined finite-size transition point. 

In the canonical ensemble it is impossible to define an order parameter which 
exhibits the typical features of spontaneous symmetry breaking already in 
finite systems. Usually, a pseudo critical temperature is defined via the position 
of the maximum of some thermal quantity. However, this definition introduces 
an ambiguity as different quantities, e.g. the specific heat or the susceptibility, 
normally have their maxima located at different temperatures. 

In order to determine relevant microcanonical quantities with high precision, 
very accurate estimations of the density of states are needed. The data pre- 
sented here have been obtained with a very efficient algorithm [6] based on 
the concept of transition observables [10]. From these observables derivatives 
of the entropy with respect to the energy and/or magnetization are easily 
computed. This is a point of utmost importance as it is not the entropy itself 
which enters into the microcanonical analysis, but its derivatives with respect 
to E and M. 



The microcanonical finite-size scaling theory proposed and discussed in the 
following is formulated in such a way as to take advantage of the existence of 



a well-defined transition point in finite microcanonical systems. This leads to 
the following scaling ansatz 



L^^'""' msp,N i£c,N 



W 



(Cie,,^-e)L'/^^) 



(5) 



for the microcanonically defined order parameter in the limit of small scaling 
variables x = {Sc^n — e) L^/'^^ . C is a non-universal, model-dependent metric 
constant and H^ is a universal scaling function for a given universality class. 
Note that the microcanonical critical exponents in equation (5), which describe 
the behaviour of various quantities with respect to the specific energy, are not 
identical to the usual canonical exponents. For the order parameter one has 
j3e = (3 /{I— a), whereas the microcanonical correlation length critical exponent 
is given by z/^ = z^/(l — a) [11,5]. Here f3, a, and u are the canonical critical 
exponents. As the microcanonical order parameter varies like a square root in 
the vicinity of e^A^ for all finite system sizes A^ [5,7] the scaling function W is 
asymptotically given as a square root W{x) ~ ^/x for small scaling variables 
X. One remarkable feature of equation (5) is the absence of any non-universal 
quantity related to the infinite system. Especially, the location of the infinite 
system critical point £c,oo does not enter in the definition of the scaling variable 

X. 
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Fig. 2. Microcanonical finite-size scaling plot for the three-dimensional Ising model. 
The values (3e/i^s = 0.54 ± 0.03 and l/u^ = 1.43 ± 0.04 result from the best data 
collapse. 

In Figure 2 we test the ansatz (5) by plotting L'^^/'^^ '"^sp,Ar as a function of 
X = {£c,N — ^) L^^'^^ for several, altogether rather small, three-dimensional 
Ising models. Here, L ranges from 4 to 16. Using a recently proposed method 
for quantifying the nature of a data collapse [12], the optimal exponents and 



their error bars can be obtained. For the 3d Ising model our small system data 
yield the values (ie/ve = 0.54 ± 0.03 and l/z/^ = 1.43 ± 0.04, in remarkable 
agreement with the expected values P^/ve = 0.52 and l/u^ = 1.43. For the Gl 
model, a similar study for system sizes ranging from L = A to L = 12 yields 
the values jSe/ui; = 0.51 ±0.03 and l/ue = 1.43 ±0.04, which again agrees with 
the literature values. For the 2d three state Potts models system sizes ranging 
from L = 6 to 18 have been investigated. The best data collapse results in the 
estimates l/z/^ = 0.79±0.03 and Pe/iyg = 0.13±0.03 for the critical exponents. 
The agreement with the exactly known values l/ue = 4/5 and Pe/^e = 2/15 
is again very good considering the smallness of the systems. 
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Fig. 3. Micro canonical finite-size scaling plot both for the Ising (I) and for the 
Gl model. Adjusting the non-universal metric constants, see equation (5), with 
CqiJCi = 0.219, the data points of both models fall on a unique master curve, 
thus demonstrating the universality of the finite-size scaling function W . The inset 
shows the scaHng functions for the Ising model (solid curve) and the 2d Potts model 
(dashed) . The scaling variable x of the Potts model is rescaled so that the amplitude 
of the square root function in the limit x — > are identical (compare relation (5)). 
The curves are only displayed for the range of the scaling variable for which data 
for the Potts model have been obtained. 

In Figure 3 we investigate whether the scaling function W is indeed universal. 
Adjusting the values of the non-universal constants C/ and Cqi-, the data of 
both models should fall on a common master curve for models within the 
same universality class. With Cgi/Ci = 0.219 a unique curve is observed, 
thus demonstrating the universality of the finite-size scaling function W. In 
the inset of Figure 3 the scaling functions of the Ising model and of the Potts 
model are compared. As the two model systems belong to different universality 
classes different scaling functions are expected, which is indeed confirmed by 
the data. 



Our approach differs from earlier attempts at a microcanonical finite-size scal- 
ing theory [5,13,14,15,16] in various regards. In Ref. [5,14] it was supposed that 
the entropy of finite systems was a homogeneous function in the vicinity of 
the transition point £2,00 of the infinite system. The resulting scaling relations 
involve ^coo in the scaling variable and lead to rather poor results for the range 
of system sizes considered in the present work. However, the scaling relations 
deduced in [5] are expected to be valid in the asymptotic regime L —>■ 00. 
The authors of Ref. [16] studied the enthalpy instead of the entropy and de- 
veloped a finite-size scaling theory in complete analogy to the canonical case. 
In the same way as in the canonical ensemble they defined a finite-size order 
parameter which differs from zero for all energies. No noticeable differences 
with canonical results were found in that approach. 

In conclusion, we have presented numerical evidence that a suitable micro- 
canonical finite-size scaling ansatz with the correct values of the critical ex- 
ponents leads to good data collapse of microcanonical finite-size quantities 
of moderately small systems. It follows from this microcanonical ansatz that 
critical exponents can easily be obtained from finite microcanonical systems 
without any a priori knowledge of infinite system quantities. The expected 
universality of the scaling function for a given universality class has been 
demonstrated. The microcanonical finite-size scaling theory, presented for the 
order parameter in the present work, can also be formulated more generally 
in terms of scaling relations of the entropy of finite systems considered as a 
function of the energy, the magnetization and the inverse system size [17]. 
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